Superconducting proximity effect in quantum wires without time-reversal symmetry 
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We study the superconducting proximity effect in a quantum wire with broken time-reversal (TR) 
symmetry connected to a conventional superconductor. We consider the situation of a strong TR- 
symmetry breaking, so that Cooper pairs entering the wire from the superconductor are immediately 
destroyed. Nevertheless, some traces of the proximity effect survive: for example, the local elec- 
tronic density of states (LDOS) is influenced by the proximity to the superconductor, provided that 
localization effects are taken into account. With the help of the supersymmetric sigma model, we 
calculate the average LDOS in such a system. The LDOS in the wire is strongly modified close to 
the interface with the superconductor at energies near the Fermi level. The relevant distances from 
the interface are of the order of the localization length, and the size of the energy window around 
the Fermi level is of the order of the mean level spacing at the localization length. Remarkably, the 
sign of the effect is sensitive to the way the TR symmetry is broken: In the spin-symmetric case 
(orbital magnetic field), the LDOS is depleted near the Fermi energy, whereas for the broken spin 
symmetry (magnetic impurities), the LDOS at the Fermi energy is enhanced. 

PACS numbers: 74.45.+C, 73.20.Fz, 73.21.Hb 



I. INTRODUCTION 

Proximity effect in normal-superconducting structures 
is a phenomenon of induced superconducting correlations 
in the normal metal (N) in an electric contact with a 
superconductor (S)ii Such correlations arise due to An- 
dreev reflections at the NS interface: an electron from 
the N part reflects as a hole by emitting a Cooper pair 
into the S part£. The extent of such proximity correla- 
tions is determined by the structure and geometry of the 
N part: in a small N grain, the correlations are uniform 
over its volume, while in a large N contact the proximity 
correlations extend over some distance determined by the 
energy of the electrons relative to the Fermi energyj 3 ^ 

One of the signatures of the proximity effect is the 
modification of the local electronic density of states 
(LDOS) in the N part. Such a modification is most pro- 
nounced in the case of chaotic electron dynamics in the N 
part, when a so-called "minigap" is formed in the density 
of states £ The size of the minigap can be estimated as 
min(A, l/r esc ), where A is the superconducting gap and 
r csc is the time required for an electron in the N region to 
establish a contact with the superconductor— (hereafter 
we assume h = 1). In particular, for a long diffusive wire 
(with its length L exceeding the superconducting coher- 
ence length) and transparent NS interface, the minigap 
is of the order of the Thouless energy,— 

E Th = D/L 2 , (1) 

where D is the diffusion coefficient. 

The proximity effect is sensitive to the time-reversal 
(TR) symmetry in the N part, which is necessary for 
superconducting correlations. If the TR symmetry is 



broken (e.g., by a strong magnetic field or by magnetic 
impurities), then the conventional quasiclassical theory 
predicts that no proximity effect can survive beyond the 
distance over which the TR symmetry is broken^ 

In this paper, we show that this picture is not com- 
plete: the proximity effect does survive even in the ab- 
sence of the TR symmetry, although in a weaker form, 
if one takes into account Anderson localization. Namely, 
in the localized phase with the localization length £, we 
expect that the average LDOS is modified on the length 
scale £ away from the NS interface in the energy window 
of the order 

^ = D/e (2) 

around the Fermi energy. In the limit of a small N grain, 
this LDOS modification crosses over to the well-known 
result for the density of states in the random-matrix the- 
ory (RMT) in superconducting symmetry classes^ in this 
limit, the energy window where the density of states is 
modified becomes of the order of the energy-level spacing 
in the grain. 

While we expect that this qualitative picture is valid 
for a quite general geometry of the NS system, in the 
present paper we support it with an explicit calcula- 
tion for the quasi-one-dimensional geometry (Fig. [TJ. 
Namely, we consider a quasi-one-dimensional wire (of a 
finite length L) with a broken TR symmetry in a contact 
with a superconductor. Using the method of the non- 
linear supersymmetric sigma- model ; 10 ! 11 in conjunction 
with recent exact results for localization in quasi-one- 
dimensional unitary wires,— ~— we calculate the average 
LDOS as a function of energy and coordinate along the 
wire in different limiting cases. 
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superconducting correlations decay into the normal re- 
gion at the diffusive length scale 



FIG. 1: (Color online) Quantum wire of length L coupled to 
a massive superconductor at x — 0. We calculate the average 
local density of states {pe(x)) for an arbitrary relation be- 
tween x, L, and the localization length £. The dimensionless 
distances t and s are used in Eq. (|13|) . 



A remarkable detail of our analysis is the two different 
ways of breaking the TR symmetry. Namely, it may be 
broken either without breaking the spin-rotational sym- 
metry (e.g., by a strong magnetic field inducing a TR- 
symmetry-breaking vector potential, but only a negligi- 
ble Zeeman field) or with breaking the spin-rotational 
symmetry (e.g., by including magnetic impurities cou- 
pled to the spins of electrons). These two possibilities 
correspond to different symmetry classes in the RMT 
classification (C and D, respectively)^ and exhibit quite 
different types of the proximity effect. First, while in 
class C the LDOS is suppressed at low energies and at 
short distances, in class D it is enhanced (similarly to the 
RMT results). Second, in long wires (£ 3> £) at small 
energies (E <C A^), the proximity effect extends to the 
Mott length scaled 



L M = 2£ln(A s /E) 



(3) 



in class C, while only to the localization length £ in class 
D (Fig. [5]). Such a behavior is related to the presence 
(absence) of repulsion between an energy level and its 
mirror counterpart in the symmetry class C (D) and is 
discussed in Sec. IVII 

The paper is organized as follows. In Sec. [II] we explain 
the relation of the proximity effect in the absence of the 
TR symmetry to Anderson localization. Evaluation of 
the LDOS for a quasi-one-dimensional wire with the bro- 
ken TR symmetry is reduced to the unitary sigma-modcl 
in Sec. IIII1 and its exact solution is presented in Sec. IIVI 
Resulting expressions for the LDOS in various limiting 
cases are derived in Sec.fV] Our findings arc summarized 
in Sec. IVII Technical details are relegated to several Ap- 
pendices. 



II. PROXIMITY EFFECT WITHOUT 
TIME-REVERSAL SYMMETRY: 
ROLE OF LOCALIZATION 

In a diffusive system (without TR symmetry breaking), 
the superconducting proximity effect is described by the 
quasiclassical Usadel equation-^. For a given energy, E, 



L E = y/D/E. 



(4) 



The TR symmetry is needed to establish particle-hole 
correlations (the soft "cooperon" modest) in the wire. 
In this paper, we consider the opposite situation: the 
proximity effect in the case of the broken TR symmetry. 
It may be broken in two different ways: by an exter- 
nal magnetic field (symmetry class C) or by magnetic 
impurities (symmetry class D). In both cases, the spec- 
trum of the cooperon modes acquires a gap leading to 
their exponential decay at some characteristic scale L c . 
This scale is set either by the magnetic length or by the 
spin diffusion length. As long as L c is shorter than Le 
(at sufficiently low energies and, e.g., strong magnetic 
field), superconductive coherence brought into the nor- 
mal metal by cooperon modes exists only in a thin layer 
of length L c near the boundary, decaying exponentially 
at larger distances. 

This exponential decay of superconductive correlations 
at L 3> L c follows from the Usadel equation^ The lat- 
ter is an effective tool for nonperturbative summation of 
tree-like diagrams^ but neglects loop corrections respon- 
sible for quantum localization. We will show below that 
these corrections, though suppressed as 1/N (N is the 
number of the conducting channels in the wire) , give rise 
to the proximity effect in the absence of the TR symme- 
try. 

Localization length in a N wire with the broken TR 
symmetry is given b y 11 ' 19 ! 20 



£ = 2nvAD ~ Nl. 



(5) 



Here A is the wire cross section, I is the mean-free path, 
and v is the bulk density of states, whose definition de- 
pends on spin degeneracy: 

• In class C (an external magnetic field), the spin 
symmetry is preserved, and hence all electron states 
are doubly degenerate. In this case, the density of 
states v [in Eq. ([3]) and throughout the paper] is 
defined per spin projection. 

• In class D (magnetic impurities) , the disorder mixes 
spin components, and v is defined as the total den- 
sity of states including spin. 

Relation between the proximity effect in the absence of 
the TR symmetry and localization can be visualized dia- 
grammatically in the limit of poor transparency of the SN 
interface when all transmission coefficients of the barrier 
are small, T<1. Then the influence of the superconduc- 
tor can be treated perturbatively. In the TR-symmetric 
case, the leading process is shown in Fig.[2Ja): A Cooper 
pair tunnels at a point ri, propagates to the observation 
point r, and returns back to the superconductor at a 
point T2- Integrations over r^a taken across the wire sec- 
tion are independent, and the resulting correction to the 
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We define the quasiparticle LDOS normalized to the 
bulk value as 



FIG. 2: Diagrams for the proximity-induced correction to 
the local density of states (pe(y)) to the lowest-order in the 
tunneling transparency of the SN interface: (a) TR-symmetric 
case, when Cooper pairs independently tunnel at ri and T2 
and propagate in the N part as cooperons (shadowed); (b) 
broken TR-symmetry case, when Cooper pairs tunnel at ri ~ 
r2, propagating further as two diffusons (shadowed). 



LDOS (5pe{y)) is proportional to N 2 T 2 . In the case of 
broken TR symmetry, cooperons are suppressed, but the 
LDOS is still affected by the process depicted in Fig.[2Jb): 
A tunneling Cooper pair needs to be converted to a pair 
of diffusons which can reach the observation point. For 
such a process, the tunneling coordinates nearly coincide, 
ri ~ r2, and integration over them brings the first power 
of the wire cross section: (<5ps(r)) oc NT 2 . At the same 
time we see that the spatial behavior of the LDOS cor- 
rections (a) and (b) are identical: both decay at the scale 
Le- The \/N oc l/£ suppression of the LDOS correction 
in the absence of the TR symmetry is reminiscent of lo- 
calization and is related to one-loop structure of the dia- 
gram (b) , contrary to the tree-like diagram (a) accounted 
for by the Usadel equation. 



III. REDUCTION TO THE UNITARY 
SIGMA-MODEL 

Electronic states in superconducting systems are de- 
scribed with the help of the Bogolyubov-de Genncs 
(BdG) Hamiltonian, which acts as a matrix in the Nambu 
space: 



HBclG 



H A 

A* -QWQ- 1 



(0) 



where % is the single-particle Hamiltonian, A is the order 
parameter field, and O is a unitary matrix that defines 
the time-reversal operation: ip \— > Qtp* . If the spin sym- 
metry is preserved (class C), we write the BdG Hamil- 
tonian for one spin projection (e.g., spin- up electrons 
and spin-down holes) and fix O = 1. On the contrary, in 
class D, when the single-particle spin dynamics is non- 
trivial, we include both spin projections and put O = is y 
(the Pauli matrix in the spin space). Thus, in class D the 
dimension of the Hamiltonian ([5]) is twice larger than in 
class C. 



^(r) = -^Imtrgf(r,r) 



(7) 



It is expressed in terms of the retarded Green function of 
the BdG Hamiltonian (O, 



(8) 



Note that the definition of v in Eq. (JT]) is different in the 
C and D classes (see Section [H]), which is consistent with 
the difference in the dimensionality of T-Lbag- 

In a field-theoretical language^ - — disorder averaging 
of pe{y) is performed by representing as a super- 
symmetric functional integral and following the standard 
line of the sigma-model derivation^. The details of the 
derivation are presented in Appendix |Al and below we 
summarize the results. 

In the normal part (A = 0) of a hybrid NS system, the 
Nambu-Gor'kov Green function G§ essentially involves a 
pair of the retarded and advanced normal-metal Green 
functions with opposite energies, G E and G A : E , 

which get coupled due to Andreev reflection off the 
superconducting order parameter. Therefore the sigma- 
model for (Ge) m the normal region of an NS system can 
be exactly rewritten in terms of Evetov's sigma-model for 
the product (G^G^p). This relation has been recently 
demonstrated in Ref. l26l where an explicit mapping be- 
tween the two models was constructed. In our problem, 
(pe(^)) may be calculated in terms of the usual Efetov's 
sigma model written for the normal wire, but supple- 
mented with a boundary condition at the NS interface 
responsible for Andreev reflections. 

The TR symmetry in the wire can be broken either by 
the orbital magnetic field or by magnetic impurities. In 
the former case each level is double degenerate due to the 
spin symmetry, whereas in the latter case this Kramers 
degeneracy is lifted. For sufficiently strong symmetry 
breaking^ (or at length scale larger than the length L c 
associated with the TR symmetry breaking) Models Ha 
and lib in Efetov's classification^ are realized. In both 
cases, the cooperon degrees of freedom are frozen out, 
and the resulting sigma model is written in terms of a 
4x4 supermatrix acting in the Fermi-Bose (FB) and 
retarded-advanced (RA) spaces. 

For a normal metal, the Models Ha and lib are math- 
ematically equivalent (unitary sigma model). This is not 
the case for the normal part of a hybrid system, since 
the form of the effective boundary condition at the NS 
interface is sensitive to the way the TR symmetry is bro- 
ken. In the presence of Andreev scattering, the Models 
Ha and lib correspond to the symmetry classes C and D 
in the classification of Ref. 

Thus, in the limit of strongly broken TR symmetry, the 
proximity effect in the normal part of an NS system is 
described by Efetov's unitary sigma-model. The average 
(normalized) LDOS is given by the functional integral 
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over the normal-metal region: 

( Pe (t)) = \ Re J str(k^Q(r))e- slQ] - SrlQ] DQ. (9) 

The action of the model is separated into the bulk diffu- 
sive part S[Q] and the boundary term Sr[Q], derived in 
Appendix [A] The bulk action has the standard form for 
the unitary class sigma model, 



s = T / dr str [• D ( VQ ^ 2 + UEA Q] 



(10) 



with the 4x4 supermatrix Q acting in the FB and RA 
spaces. The matrix A = af- A is the metallic saddle point , 
and the supersymmetr y b reaking matrix k = (we 
follow notations of Ref . fill) . As before, D is the diffusion 
coefficient and v is the density of states at the Fermi level 
per one spin projection if the spin is conserved (class C) 
and including both spin projections in the case of broken 
spin symmetry (class D). 

The boundary action is derived in Appendix [A] (see 
also Ref. |26| ). At low energies E <C A, it has the form 

S r = -i^strln[l - e^Q(0)EQ T (Q)E T }. (11) 



Here the parameters are related to the transmission 
coefficients at the SN interface: Tj = l/cosh 2 /3i, with 
i labeling open channels. In the case of conserved spin 
(class C), the channels are double degenerate but the in- 
dex i counts them only once. With broken spin symme- 
try (class D), this degeneracy is lifted. In both cases, the 
normal-state conductance can be written as a product 
GogN, where c/jv = J2i Ti, and the conductance quantum 
is defined as Go = 2e 2 /h for class C and Go = e 2 /h for 
class D. 

The main ingredient which distinguishes between the 
symmetry classes C and D is the matrix S entering 
Eq. (fTTj) . For class C it has been derived in Ref. H^, and 
derivation for both classes is presented in Appendix IA1 



class C: 



class D: 



_RA 

V 





_RA 









T RA 

J x 



_RA 



FB 



FB 



(12a) 
(12b) 



The matrix H has a nontrivial structure in the super- 
space, which is the mathematical reason why Eq. © re- 
sults in a nontrivial LDOS. 

Note that the matrix S does not belong to the stan- 
dard unitary manifold for Q matrices: it does not obey 
the condition S 2 = 1 [in the F (B) sector for class C (D), 
respectively] . One finds that the saddle-point solution for 
the action 5[Q]+5r[Q] is simply Q(r) = A, leading to the 
metallic density of states, p(r) = 1, indicating no prox- 
imity effect at the level of the Usadcl equation. However 
proximity effect absent in the quasiclassical approxima- 
tion will manifest itself once fluctuations around A are 



taken into account. The same fluctuations are responsi- 
ble for localization. Thus the superconducting proximity 
effect in a normal metal with broken TR symmetry is 
inevitably related to localization, leading to strong mod- 
ification of the LDOS near the SN interface at energies 
E < A^ and length scales x < £. 



IV. EXACT SOLUTION FOR QUASI-1D WIRES 

Here we apply the general framework described in the 
previous Section in order to calculate the average LDOS, 
(pe(x)), in a finite qlD wire of length L coupled to a 
superconductor at x = 0, see Fig. [TJ We will work in the 
subgap limit, E <C A, and assume that the TR symmetry 
in the wire is completely broken. 

The one-dimensional sigma model can be solved 
exactly by mapping onto effective quantum mechanics,— 
with the x coordinate playing the role of the imaginary 
time. In this formalism, evaluation of the functional in- 
tegral © is reduced to solving Schrodinger equations 
for the wave functions ^(Q) on the sigma- model man- 
ifold. This technical procedure described in Appendix IA1 
leaves us with the object ^>(Xp,Xb) depending only on 
the "eigenvalues" Xp and As of the Q matrix. 

The main ingredient for calculation of (pe{x)} is the 
wave function ty(\p, \s',t, s) which accounts for coher- 
ent motion of the electron and hole in the wire. It can be 
obtained by successive application of two evolution oper- 
ators on unity (corresponds to open boundary conditions 
at the free end of the wire) : 



V(\ F ,\ B ;t,s) = e 



-2Ht-2Hs 



ol, 



where (see Fig. [T]) 

t = x/S, s = (L- x)/t 



(13) 



(14) 



The Hamiltonians H and H govern imaginary-time evo- 
lution at the segments [0, x] and [a;, L], respectively. They 
have the form2£ 



H 



(Xb-X f ) 2 



d 



1-A| 



d 



dX F (A s - X F ) 2 dX F 



d X 2 B -l 







— (X B -X F ), (15) 



dXs (Xb - X F ) 2 8Xb 
andi^ 

H = (Xb — Xp) 1 H(Xb — Xp) = Hb + Hp 

with 



Hp 



1 d . 2 d k 2 

2 9Al ( " + IT"' 

1 a -(l-A^> 9 



2dXi 



^X 
dX F 16 F ' 



(16) 

(17a) 
(17b) 
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where the dimcnsionless quantity k stands for 

8iE 
'AT 



,.2 



(18) 



[it is consistent with the notations of Refs. [l2| and [Tj, 
with the frequency w = E — (— E) = 2E]. 

Unfortunately, the function ^(A^, Xb] t, s) cannot be 
generally obtained in a closed form. The only exception 
is the zero mode of the Hamiltonian H, ^o(^f,Xb) = 
^(Xf, Xb] 0, oo), given byH 

*o(A F , X b ) = h{q)pKi(p) + qh(q)K (p), (19) 

where 

p = kV(A b + 1)/2, q = Ky/{\ F + \)/2. (20) 

The function ^(Xf, Xb', t, s) should be finally inte- 
grated over Q(0) with the weight e~ Sr ^ in order to ob- 
tain (pe{x)}. At this stage calculations for the symmetry 
classes C and D are different due to a different form of 
the superconducting matrix f|12[>. 



A. General expression for class C 

For the symmetry class C, calculations presented in 
Appendix [B] lead to the following exact expression: 



(p E (x)) =1 + 1; Re J dX F J dX B 



de~ s ° ( Xb ^ 

x^— *(A F ,A B ;t,8). (21) 

oXb 

Parameters of the SN interface are encoded in the bound- 
ary action: 

S c (A B ) = ^ln[l + 7I(A|-l)], (22) 

i 

where % = T 2 /(2 — Tj) 2 is the Andreev transmission of 
the z'th channel^ 

The boundary action (|22|) suppresses fluctuation of the 
bosonic variable A_b. The strength of the coupling to the 
superconductor is characterized by the dimcnsionless [in 
units of Go defined below Eq. (jlip] Andreev conductance 
of the interface, 



9 A 



(23) 



In the limit of gA S> 1, the integral over As in Eq. (|2lj) 
comes from As — 1 < 1/gA- If ^(Xf, Xb; t, s) is a slow 
function of As in this region then Eq. (|2"Tj) simplifies to 
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FIG. 3: (Color online) Regions of different behavior of the 
average LDOS in the coordinates L/f, E/A$. 



B. General expression for class D 

In close analogy with Eq. (|2"Tj) . calculations in Ap- 
pendix [B] lead for the symmetry class D: 

"1 p oo 



(p E (x)) = l + ^R c J dX F J dX 



dXi 



*{\ F ,\ B ;t,s), (26) 



with the boundary action 

S$(\ F ) = -^U[l-Ti(l-\ F )]. (27) 

i 

For large gA [Eq. (f2~5j) ]. Xf is pinned to ±1 and we obtain 
1 f°° 

(p B (x)) = 1 + -Re / dX B [^(l,XB]t,s) 

-*(-l,A B ;t,«)]. (28) 

V. RESULTS 

Here we analyze the general expressions (|2T|) and ([26]) 
in various regions of the system parameters, schemati- 
cally shown in Fig. [3] We will be mainly interested in 
the limit gA — > oo [Eq. ([25]) ] when the LDOS is given by 
Eqs. (|M|) and (f2"5)l , and present the result for an arbitrary 
gA only in the perturbative regime discussed in Sec. IV Bl 



1 r 1 

(p E {x)) = 1 - - Re / dX F ^(X F , 1; t, s). (24) A. Random-matrix regime (L < £ and E < £ T h) 

2 7-i 



As we will see below, approximation (f2"4"|) applies if 

g A »mm(jZ/L Et E/5,i) t (25) 
where S = (2vAL)~ 1 is the mean level spacing in the 



We start the analysis of the general expression (f2"4"| for 
the LDOS with the simplest case of short wires, L -C £. 
In the limit of _E <C £rh, electrons have enough time to 
explore the whole available space and one should recover 
the RMT statistics. In this case, one can neglect the 



6 




FIG. 4: (Color online) Energy dependence of the RMT den- 
sity of states for the symmetry classes C and D. 



terms with derivatives in the Hamiltonians: H ss H 
(k 2 /16)(Ab — Xf), and one simply gets 



V(X F ,X B ;t,s) = e -(K 2 /8)(A B -A F )(L/0 



(29) 



independently of x. Substituting into Eqs. (|24|) and (|28|) 
one readily recovers the RMT results for the symmetry 
classes C (minus sign) and D (plus sign)£ 



(pe(x)) = 1 T 



sin(27r£y<5) 



(30) 



with the mean level spacing 6 = (2vAL)~ 1 [this definition 
of the level spacing takes into account both electron and 
hole states, hence it contains an additional factor 1/2]. 
Energy dependence of the RMT density of states is shown 
in Fig. |H 

The gradient terms in the sigma model may be taken 
into account perturbatively, which results in a position- 
dependent correction to the RMT density of states (|3U1) : 



(Sp B (x)) = ± 



L 



it 



sm{2irE/5) 



2-rrE/S 

O(L/0 + O(E/E Th )\, (31) 



1 t- cos(2tt£/<5) 



with the upper (lower) signs corresponding to the sym- 
metry class C (D). Note that in the symmetry class C 
both the main result (|3"0"]) and the correction (|3"Tj) vanish 
at zero energy. The gradient terms in the Hamiltonians 
H and H become significant if the length of the wire L 
exceeds £ or at high energies E > £rh- This establishes 
the boundaries of the RMT regime, see Fig. [3J 



B. Perturbative regime (£>ETh,Aj) 

For E > max^Th, Ag), the LDOS is close to 1 ev- 
erywhere in the wire, with the difference (pe{x)) — 1 de- 
caying at the length scale L B from the SN boundary. In 
this case only small deviations of Xf and Xb from 1 arc 



important. Then the Hamiltonians H and H can be sim- 
plified in the vicinity ofA F = As = l ("north pole" of 
the fermionic sphere) yielding an effective linear oscilla- 
tor model. In this approximation, the wave function (|13[) 
takes the following form: 



*(A F ,A s ;i, 



x cxp 



cos1i 2 (ks/2) 
cosh 2 («:(t + s)/2) 

-" (AB ' Af) tanh K(t + - s) 
4 2 



(32) 



The density of states for the two symmetry classes then 
follows from Eqs. (|2"Tj) and (j2l))) . where only small de- 
viations of Xf and Xb from 1 arc important. In this 
limit, the boundary actions (f22|) and (|2T[) take the form 
S§(X B ) = ( 5 A/2)(A B -l)andS* D (A F ) = (g A /2)(l- X F ). 
Integrating over 1 — Xp and Xb — 1, we obtain 



( PE (x)) = 1 =f 4 Re 



cosh 2 ks/2 
KsinliK(£ + s) 
9a 



g A + (k/2) tanh(«;(t + s)/2) ' 



(33) 



As in Eq. (|30p . the minus (plus) sign corresponds to the 
symmetry class C (D). 

Modification of the metallic LDOS by the proximity ef- 
fect is maximal for good SN contacts with g A gN(LE), 
where <7at(£_e) = ^/L B = WE/A^ is the normal-state 
conductance of the wire of length Le- In this case, the 
second factor in the correction (|33p may be approximated 
by one, and the result (in the physical units) reads 

{pe(x)) = l T Re ™* a jq -g^ - , (34) 



(l-i)amh[2(l-i)L/L B ] 



In the particular case of an infinite wire, Eq. (|34[) simpli- 
fies to 



(Pe{x)) 



Le o„it ( 2x 7r 
1 T -A e- 2x ' LE cos — + - 
V2£ \L E 4 



(35) 



[Similar simplifications are also possible in the opposite 
limit g A <C gi\f(LE): in that case, the perturbative cor- 
rections in Eqs. (|34p and (|35p get multiplied by g A and 
correspond to the diagram in Fig. [2Jd.] This result agrees 
with our argument in Sec.|Hj the proximity effect extends 
to the normal region at the length scale Le (just like in 
the conventional Usadel equation), but its magnitude is 
small as L B /£ *C 1 (in the perturbative regime) . 

Note that our perturbative calculation gives only the 
contribution to the LDOS from the vicinity of the north 
pole of the fermionic sphere (Xp = Xb = 1). In 
fact, the expressions ([3"2"j) - (fM]) provide a valid pertur- 
bative treatment of the north-pole contribution at en- 
ergies E 3> max(<5, Ag). In particular, in the window 
6 <C E <C Et\i-, the north-pole contribution repro- 
duces the non-oscillating parts of the RMT result (|30p . 



E/A e 



FIG. 5: (Color online) Energy dependence of the LDOS in the 
vicinity of the NS boundary, {p^ (0)), for the semi-infinite 
wire. The lower (upper) curve corresponds to the symmetry 
class C (D). 



The south pole of the fermionic sphere (Xp = —1, 
Xb = 1) is an alternative saddle point of the action, 
which is responsible for the terms oscillating in energy 
with the period S in the RMT [see, e.g., Eqs. ([30|) and 
(|3~Tj) ].— At energies E ^> i?Th, its contribution is ex- 
ponentially suppressed, and therefore it was neglected 
in our perturbative calculation. However, at energies 
E < -Et^ the south-pole contribution exceeds, by abso- 
lute value, the correction to unity in Eq. Note that 
at E 3> max(<5, A^) the south-pole contribution may also 
be found in a perturbative scheme analogous to the one 
above, but expanding around the south pole instead. We 
do not report this calculation here, but only remark that, 
within the window 8 <C E <C -Etiu it can reproduce the 
oscillating terms in the RMT expansion (f30| and (f3~T|) . 
In other words, this energy window represents an over- 
lap between the RMT and perturbative approaches (if 
the south pole is taken into account). 



C. Localized regime (L>^ and E <C A^) 

Finally, we turn to the limit of long wires, L>^. Here 
the difference between the symmetry classes C and D 
becomes more pronounced than just the sign in Eqs. (|30j) 
and ([33]) . and they require separate considerations. As 
in subsection A, our results in this section are restricted 
to the "good contact" limit gA > 1 [sec Eq. (|25l) ]. 



1. Symmetry class C 

Using Eq. (|19p for the zero mode, one immediately 
obtains an exact expression for the energy dependence of 
the LDOS at the SN interface for the semi-infinite wire: 

( PE (0)) = I -2 Re [Ii («)#!(«) + I 2 (k)K (k)] , (36) 



where the imaginary parameter n 2 is defined in Eq. (|18[) . 
The plot of this function is shown in Fig. [S] (lower curve). 



FIG. 



6: (Color online) Dependence of the low-energy (E <C 
15; LDOS, (pg\x)), on the distance from the SN interface 
for the semi-infinite wire. In the symmetry class C (lower 
curve, Eq. (|39[1 . sketch), LDOS depletion extends to half of 
the Mott scale (O, while in the symmetry class D (upper 
curve, Eq. (|42[) k excess LDOS accumulated near the interface 
relaxes at the localization length. 



Its asymptotics are 



(Mo)) 



'ttE/2A^ 



1 - y/A^/AE, £>A S . 



(37) 



[The E ^> A^ asymptotics agrees with Eq. (|35|) obtained 
by perturbative methods.] 

For small x <C £, the function (pe(x)) can be obtained 
as a power series in x by expanding the evolution oper- 
ator e" 2i?t in t = x/£ (cf. Sec. VI of Ref. [3). In the 
nonperturbative limit of E <C A^ , this procedure yields 



(Pe(x)) 



nE 
2A7 



1 



2^ + 2^ 



(38) 



According to Eq. (|3"8"|) . for small energies, E <§; Aj, the 
LDOS is strongly suppressed even at the distances from 
the SN interface comparable to the localization length, 
x ~ £. In this nonperturbative regime, the LDOS deple- 
tion propagates to a much larger distance of the order of 
the Mott scale ([3]). Indeed, in the limit E <§; A^, the wave 
function ^(Q) is uniformly spread over the Fermionic 
sphere, and Eq. ([M]) yields (pe(x)) = 1 - *(1, 1; t, s). To 
the leading approximation, 'F at the north pole is calcu- 
lated in Appendix [Cl and one gets (for \t— iM /2| <C £a/): 



(pe(x)) w - + - erf f 



2y/t J 



(39) 



where tu = Lm/£- This dependence is shown in Fig. [6] 
(lower curve). 



2. Symmetry class D 

Substituting the zero mode (fT!?|) into Eq. (f?5|) . one 
readily obtains the energy dependence of the LDOS at 
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the SN interface for the semi-infinite wire: 

(pjs(O)) = 1 + 2Rd[I 1 (k)K 1 (k) + (J (k) - 1) K 2 (k)} , 

(40) 

with the asymptotic behavior 



(pe(0)) 



'3-ttE/A ( , £«A c , 



1 + s/Aj4E, E » A € . 



(41) 



[As in class C, the E ^> asymptotics agrees with the 
perturbative result (|35l) .] 

The function (pe(0)) is shown in Fig. [5] (upper curve). 
Remarkably, at the lowest energies, E <C A^, the LDOS 
at the SN boundary is enhanced by the factor of three 
compared to the normal case. This behavior should be 
contrasted with the LDOS behavior in class C, where the 
LDOS near the SN interface is depleted at low energies. 

At zero energy, E = 0, the LDOS can be obtained ana- 
lytically for an arbitrary relation between x, L and £. In 
the limit k — > 0, the Hamiltonians H and H get simpli- 
fied, which allows us to evaluate the operator exponents 
in Eq. (|T5|) with the help of the Lebedev-Kontorovich 
transformation 3 ^, see Appendix |Dj The general depen- 
dence of (po(x)) on the two parameters t and s given by 
Eq. (|D12[) is quite complicated, and we present here its 
limits at t = and at s = oo. 

For the semi-infinite wire, the zero-energy LDOS is 
given by 



7T d 

Idt 

o 

k dk 



-2I\ 



sinh(7rfc/2) _, 1 



(l + k z )t/4 



coslr(7rfc/2) 



(42) 



where t = x/£. The function pg°°' (x) is shown in Fig. |6] 
Near the SN interface, the LDOS exceeds the unper- 
turbed value by the factor of three. 

The excess LDOS accumulated near the boundary with 
the superconductor is relaxed to the metallic value at 
the scale of the localization length [whereas the LDOS 
depletion in the C-class case is restored at a larger Mott 
length scale, see Eq. ([39])]. 

For a finite wire of length L, the zero-energy LDOS at 
the SN interface is given by 



U i) (0))=3 + i 



°° irk 
k dk tanh — 
2 



rz-2 



7 



k 2 + 1 



-(fc 2 + l)s/4 



k 2 + 1 

fc 2 + 9 ( 



-(fc 2 +9)s/4 



(43) 



where s = This function shown in Fig. [7] interpo- 

lates between the RMT value of 2 for short wires (L <C £) 
and the value of 3 in the localized regime (i > {). For 
short wires, (p ( L) {0)) = 2 + 2L/3£ + . . . (dashed line in 
Fig- EJ) agrees with the perturbative result ([3~T]) . 




FIG. 7: (Color online) The zero-energy LDOS at the SN in- 
terface for a finite wire of length L in the symmetry class D. 
The dashed line shows the linear behavior for small L/t;: 
<Po L) (0)> = 2 + 2L/3£ + . . . , see Eq. dUJ). 



VI. DISCUSSION AND CONCLUSIONS 

In this paper, we have solved the problem of the LDOS 
in a one-dimensional wire without TR symmetry con- 
nected to a massive superconductor. Our calculation 
confirms qualitative expectations that the LDOS in the 
wire is modified in the energy window of Aj around the 
Fermi level within the distance of the order of £ (up to 
logarithmic corrections) from the interface with the su- 
perconductor. 

Propagation of superconductive correlations in the N 
region with the broken TR symmetry has been studied 
previously in the context of the Josephson coupling in a 
strong magnetic fiekh 31 ' 32 It has been shown that while 
the average Josephson coupling between a pair of S is- 
lands through the N region without the TR symmetry de- 
cays exponentially, its variance decays only algebraically. 
This phenomenon is similar to the universal conductance 
fluctuations , 33 ' 34 and originates from the diffuson-only 
diagram which are insensitive to the TR symmetry break- 
ing. In the present paper we demonstrate that proximity 
effect without the TR symmetry manifests itself already 
in the average properties, provided that localization is 
taken into account. 

The observed modification of the LDOS, in the case 
of a wire much longer than the localization length £, is 
totally due to localization effects and is absent in the 
quasiclassical description (see discussion in SectionQl]). In 
other words, localization enhances proximity effect near 
the NS interface. Qualitatively, the effect of localization 
is similar to cutting the wire at the distance £ from the 
interface, so that the density of states near the interface 
becomes similar to the RMT result in a grain of size £. 
While this similarity is only qualitative in long wires, 
we indeed find a true quantitative crossover to the RMT 
regime in the limit of short wires (much shorter than 
£). Note, however, that both in the long- wire and short- 
wire (RMT) limits, the number of energy levels involved 
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in the modification of the LDOS is of order one (since 
uA^A ^ 1). Thus the physics of this modification may 
be understood in terms of the repulsion of the lowest 
levels from their own mirror images (generated by the 
Andreev reflection), in the spirit of the superconducting 
ensembles in the RMT.— 

This picture is in agreement with another important 
observation: the sign of the correction to the LDOS de- 
pends on the symmetry of the TR-breaking terms. There 
are two possible ways to break the TR symmetry: cither 
conserving the spin symmetry (class C realizable, e.g., 
by an orbital magnetic field) or breaking it (class D re- 
alizable, e.g., by magnetic impurities). It is known from 
the RMT that corrections to the LDOS at low energies 
have opposite signs in these two symmetry classes (see 
Fig- HJ) - In class C, the levels repel from their own mirror 
images,— and therefore there is an effective repulsion from 
the Fermi level leading to a suppression of the LDOS 
around it. On the other hand, in class D, there is no 
repulsion of any level from its mirror image,— and there- 
fore the LDOS is enhanced around the Fermi level (levels 
are "squeezed" towards it). According to the results of 
our calculations, the same tendency persists also for long 
wires, away from the RMT limit (see Figs. [S] and [5]). The 
above argument about the repulsion of levels from their 
mirror images suggests that our results may be explained 
in terms of Mott hybridization 1 ^ of low-lying energy lev- 
els with their Andreev-reflected images, in the spirit of 
our recent work.— We leave this interesting question for 
future study. 

Experimentally a modification of the LDOS by one 
level is usually difficult to observe (in particular, it is 
typically not self-averaging). However similar effects may 
be of relevance in more experimentally accessible geome- 
tries (e.g., in planar NS junctions, where the corrections 
to the LDOS would be self-averaging due to a large area 
of the NS interface). Also, our calculation for the unitary 
symmetry class may provide helpful intuition for treating 
systems without TR symmetry breaking. 

An important step in studying nonperturbative local- 
ization effects in such proximity systems without TR 
symmetry breaking was done in Ref. l36l There a 
scattering-matrix formalism was used to find the mod- 
ification of the LDOS in TR-symmetric quantum wires 
close to the NS interface, in comparison with the pre- 
dictions from the Usadel equations. It was found that 
the LDOS is reduced within the energy window of the 
order around the Fermi level, which is qualitatively 
similar to our results for wires with broken TR symme- 
try. It would be interesting to extend the results of Ref. 
[3a to finite distances from the NS interface, possibly us- 
ing indications from our present work. Progress in this 
direction may also be relevant for a theoretical analy- 
sis of experiments on Josephson junctions via Anderson 
insulators i2£ 

Recently a great interest has been attracted to one- 
dimensional hybrid structures hosting Majorana fcrmions 
at the NS boundary. Such bound zero-energy eigenstates 



appear if the superconductor is of a topological type i 38 i 39 
In the random-matrix context, such systems are assigned 
to class B.— ~— An evidence in favor of Majorana exci- 
tations in an InSb nanowire coupled to a superconduc- 
tor has been reported in very recent experiments! 43 ' 44 
The zero-energy mode at the NS boundary is topolog- 
ically protected and hence robust with respect to dis- 
order. It is this robustness that makes the Majorana 
fermions promising candidates for realization of fault- 
tolerant quantum computations 4 ^ and motivates further 
intensive research in the field. 

Disorder effects on the spectral and transport char- 
acteristics of onc-dimcnsional hybrid structures of class 
B models were studied both analytically and numeri- 
cally within the scattering matrix formalism.— ~— The 
sigma-model approach developed in our paper can also 
be applied to the hybrid wires of class B. The character- 
istic distinction between classes D and B in the random- 
matrix approach is that the two disconnected parts of the 
sigma-model manifold contribute with either the same 
(class D) or opposite (class B) signs to the partition 
function.— In our onc-dimcnsional problem, this corre- 
sponds to the sign of the "south-pole" (Ap = —1) contri- 
bution in Eq. ([25)1 . By changing the difference of the two 
wave functions in Eq. (|28[) into their sum, we can readily 
infer the LDOS in class B including the spatial profile of 
the Majorana mode and the depletion of the local density 
of other low-energy states. This calculation will be the 
subject of a separate forthcoming publication* 4 ^ 
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Appendix A: Sigma model for the NS junction in 
classes C and D 

In this Appendix we outline the derivation of the sigma 
model for a normal metal - superconductor junction both 
in the case of preserved (class C) and violated (class D) 
spin symmetry. The form of the sigma-model action is 
known to be universal and depends only on the symme- 
tries of the underlying system but not on its microscopic 
details. Using this fact, we will pick a particular Hamilto- 
nian which bears all necessary symmetry properties and 
facilitates the derivation of the sigma model. 

Our strategy is as follows. We assume that the domi- 
nant disorder in the sample has the form of random po- 
tential. This potential disorder sets the elastic mean free 
path I, which is the shortest scale in the problem, and 
preserves both time-reversal and spin symmetry. This 
assumption allows us to consider both superconducting 
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and normal part of the junction on equal footing and 
describe them by a single sigma model for a diffusive 
system. Then we introduce an additional relatively weak 
source of disorder that breaks the TR symmetry in the 
normal part of the junction. This is either random vec- 
tor potential (class C) or magnetic impurities (class D). 
Violation of time-reversal the TR symmetry occurs on a 
longer scale L c 3> I. Retaining only masslcss degrees of 
freedom, we finally obtain the sigma model applicable at 
distances longer than L c . This will be the unitary class 
sigma model with boundary conditions of either class C 
or D. 

Starting point of the derivation is the Bogolyubov-de 
Genncs Hamiltonian. Its general form is given by Eq. ([6]) . 
The Hamiltonian is a 2 x 2 matrix in the Nambu space. 
We use the notation r for Pauli matrices acting in this 
space. Let us briefly discuss the structure of the Hamilto- 
nian for our particular model. The strongest disorder is 
the random potential U. Singling out the corresponding 
term, we write the Hamiltonian as 



%BdG = T-La + t z U. 
The most general form of the term Hq is 



(Al) 



class C: H = r z £(p) + r y A a • v, (A2a) 
class D: H = r z £(p) + r y A - b • s. (A2b) 

Here we have included (i) the standard kinetic term, (ii) 
the superconducting order parameter A (assumed to be 
real), and (iii) the random vector potential a (for class 
C) or the random exchange field b due to magnetic im- 
purities acting on electron spin (for class D). In the su- 
perconducting part of the junction, a = b = 0, whereas 
in the normal wire, A = 0. We proceed with the deriva- 
tion of the sigma model using the general form of the 
Bogolyubov-de Gennes Hamiltonian and will specify the 
effect of particular terms later. 

Local density of states is expressed in terms of the 
retarded Green function ([7]). This Green function can be 
generated from the path integral with the action 



S = -i J dx<P + (E + iO- HbcIg)*- 



(A3) 



Here the vector field $ contains both commuting and 
Grassmann variables and also has the structure in the 
Nambu space. 

Bogolyubov-de Genncs Hamiltonian obeys the symme- 
try H = —T y Ql-l*Q~ l Ty that provides the onset of spe- 
cific soft modes - superconducting cooperons. In order to 
include these modes into the sigma model, we rewrite the 
action (|A3[) in the doubled form, introducing the particle- 
hole structure of the fields. 

S = -if dx^[(E + iO)a z T z -T z n BdG ]^, (A4) 

1 ( $ \ . 1 



We denote Pauli matrices acting in the particle-hole 
space by a. Later, we will see that this space corre- 
sponds to the retarded-advanced structure in the Efctov 
parameterization of the unitary class sigma model. The 
supersymmetry-breaking matrix k = diag(l, — 1)fb- 

The matrix is either unity (when spin is preserved, 
class C) or is y (spin symmetry is violated, class D). In 
both cases, the vectors \1/ and are related by the iden- 
tity ^> = (C^) T with the charge-conjugation matrix C 
dependent on the symmetry class. 



class C: 



class D: 



C 



ICTy 

(7 X 



7 X 
io, 



yj fb 



(A6a) 
(A6b) 



Derivation of the sigma model proceeds with averaging 
the statistical weight e~ s over potential disorder. As- 
suming Gaussian white-noise statistics of U, averaging 
yields a quartic term ~ (^vp) 2 . This term is further de- 
coupled with the help of the auxiliary matrix field Q. 
The action acquires the form 



S = dx 



8r 



str Q 2 - [ Eu z t. 



iQ 

2r 



(A7) 

Here Ho is defined in Eqs. (|Al[) and (| A2|) . and isq is 
the density of states per one spin projection. The field 
2 is normalized such that the infinitesimal imaginary 
part i0a z r z is replaced by iQ/2r where r is the disorder- 
induced mean free time. 

Next, we integrate out the vector field and perform 
the saddle-point analysis of the resulting action 



St 



1 



S = — - / dx str Or - - str In Eo z t z + - t z Uo 



iQ. 



2r 



(A8) 

The term proportional to E as well as all the terms in %q 
other than the kinetic term [see Eq. (|A2[) ] are relatively 
small. We first consider saddle points of the action (|A8|) 
with E = and Ho = TzC(p)- There is always one sad- 
dle point equivalent to the self-consistent Born approx- 
imation for the average Green function. It corresponds 
to replacing iO in the original action with i/2r, hence 
Q = u z t z is a legitimate saddle point. We can generate 
other saddle points with the help of the global gauge sym- 
metry of the fermionic action (|A7|) . Indeed, let us change 
the variables in Eq. (|A"7|) : * i— > U^S, # H> *t7 with some 
matrix U obeying the relation UU = 1. (Charge con- 
jugation for matrices is defined as U = CU T C T .) This 
transformation cannot alter the value of the action for 
Q. At the same time, U commutes with £(p) since U is 
constant in space. Hence the rotation of $ can be com- 
pensated by Q n- UQU" 1 . This allows us to generate 
the whole saddle manifold of the sigma model, 



(A5) 



Q = Ua z T z U. 

The manifold is characterized by the conditions Q 2 
and 2=Q = CQ T C T . 



(A9) 
1 
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The sigma-model action is the result of the expansion 
of Eq. (|A8p to the first order in scalars A and E and to 
the second order in vectors VQ, a, and b. Assuming a 
and b are Gaussian white-noise random quantities, we 
obtain 

5= ir/ dxstl '{ D ( VQ ) 2 + H E °z T * + ^ A )Q 

+ L sh y (A10) 

with the symmetry-breaking term 

class C: L sb = -D(a) 2 [T z , Q] 2 , (Alia) 
class D: L sb = -^(b) 2 [r 2 s, Q} 2 . (Allb) 

The above action contains simultaneously all possible 
terms of our sigma model. Now we analyze the theory 
in particular cases keeping only the relevant terms in Eq. 
(|A10[) . Let us start with the superconducting part of the 
junction. Assuming a bulk superconductor with large 
density of states and diffusion constant, we keep only E 
and A terms and neglect fluctuations of Q. This allows 
us to fix Q = Qs at the saddle point of the sigma-model 
action, 



to b 



,11 



Qs 



( a z cos 8g sin 9g \ 
v sin6» s ~g z cos 6 S J N 



iA 

tan6> s = — . (A12) 



Now consider the normal part of the junction first with 
preserved spin symmetry (class C). Here A = while 
random magnetic field, (a) 2 ^ 0, has a strong effect vi- 
olating the TR symmetry. This leads to the additional 
constraint [r z , Q] — making the matrix Q diagonal in 
the Nambu space. The condition Q = Q with C from 
Eq. (|A6aP leads to the representation 



Q 



Q 
ZQ T Z T 



(A13) 



where 3 is the matrix defined in Eq. (|12al) . In the normal 
part of the junction, the sigma model action (|A10[) can 
be written in terms of Q with the only constraint Q 2 = 1. 
Then one obtains exactly the action (| 10[) (with v = uq) 
for the unitary class sigma model, with the particle-hole 
structure playing the role of the retarded-advanced space. 

Let us now consider the case of broken spin symmetry 
(class D). The matrix Q has a structure in the additional 
spin space and obeys the relation Q = Q with C from 
Eq. (|A6b[) . Taking into account random exchange field, 
(b) 2 ^ 0, we get the constraint [t z s, Q] = 0. It makes Q 
diagonal in the Nambu space and trivial (unit matrix) in 
the spin space. The representation (|A13[) again applies 
but now with S from Eq. (|12bp . In terms of Q, the 
sigma model action is again given by Eq. (flU)) . The trace 
over the spin space doubles the density of states in the 
prcfactor of the action: v = 2v§. 

The last ingredient of the theory is the boundary term 
pip . In terms of the matrix Q, boundary action is known 



Sr = -i^>trln(l 



QsQ)- 



(A14) 



Using Eqs. (|A12|) and (|A13|i . we trace the Nambu (and 
optionally spin) space and obtain the boundary action in 
the form 

Sr = \ str In [l + e~ 2 ^ cos 9 s {Qa z - <7 Z ZQ T Z T ) 



QZQ'E 



(A15) 



At low energies, £ < A, we have cos 0s = and obtain 
the boundary action (TlTj) . 



Appendix B: Derivation of the LDOS expressions 
(EH) and (p6j) 

1. Efetov's parametrization 

We adopt the standard Efetov's parametrizationii of 
the 4x4 matrix Q for the unitary symmetry class: 



U v 1 QoU v , 



(Bl) 



where the FB-diagonal central part, Qo, contains only 
commuting variables: 



Qo-' QF ° F 



Q BB ' ' 



/ FB 



Q FF = ' 



e-^yr^A 2 : 



-\ F 



R A 
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ie-^v/A|^l -Ai 



(B2a) 
(B2b) 
, (B2c) 



RA 



while Grassmann variables reside in the RA-diagonal ma- 
trix 



u 
Uri ~ \ v 



HA 



cr 



P U* o) FB ^ ex P Vp * 0/F 



p 



(B3a) 
(B3b) 



Integration variables change in the intervals Af £ [ — 1,1], 
Xb € [l,oo), ifiF,B € [0,27r), with the invariant measure 

d\ B d\ F dtp B dtp F (vtA s 

d ® = 71 r^2 io \2 d P d P da da - ( B4 ) 

(A b -A_f)^ (2ny 



2. Boundary term 



Substituting Eqs. (|B1[) and (|B2[) into the boundary 
action (fTTj) . after some algebra we obtain for the class C 
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[3 is given by Eq. (|12a[) 

Ab — Xf 



5 r = 



1 - 



<9A, 



5"o (As), 



(B5) 

where the action S^p (Ab) is given by Eq. ([221 . 

For class D [2 is given by Eq. (|12b[) ] we get a compli- 
mentary result: 



Sr 



X B - X F , *\r * , \ ® 
l-^^{a-p){a + P ) — 



5 D (Af), 
(B6) 



where the action S^(Xp) is given by Eq. (|27|) . 



The initial conditions for the evolution (|B9[) are Qf ~ 
g B =V(Q(x);s) at f = 0. 

Calculation of the functional integral © along the seg- 
ment [0,x] is equivalent to so solving Eqs. (|B9|) up to 
t = As a result, one ends up with a single integral 
over Q = Q(0) at the SN interface: 



(p E (x)) = Re / dQ 



sbr(kAP(Q;t,s)) 



-S r [Q] 



Using Eqs. (|B2[) and (|B8|) . we express 
str(fcAP) _ . , 



(Bll) 



(B12) 



3. Expression for the LDOS in terms of the 
effective quantum mechanics 

In the language of the effective quantum mechanics 2 ^, 
with the coordinate along the wire playing the role of 
imaginary time, calculation of the functional integral © 
for (pe(x)) is performed in two steps. First we evaluate 
the functional integral along the segment [x, L] of the 
wire, which generates the "singlet" wave function 



f(Q(i);s)=e^ s ol, 



(B7) 



where s — (L — x) /£ is the dimensionless distance from 
the open end. By symmetry, the "radial" wave function 
fy(Q(x); s) depends only on the variables X F and X B . 

At the point x, the LDOS prccxponcnt str(fcAQ(a;)) in- 
troduces the multiplet of states Q^(Q). Then, evaluation 
of the functional integral at the segment [0,x] amounts 
to following the evolution of the multiplet P = Q^{Q) 
from the observation point x down to the SN interface. 
The resulting matrix P is know n 11 ' 20 to be of the form 
(|B2[) with the elements modified according to 



Xf — y cjfXfj 
Ab — ^ 9b Ab, 



1 - K -> !f 
Ib 



1-A|, 



A B 



1 



A| 



(B8a) 
(B8b) 



The evolution equations for the functions fp and 
/b which determine the density-density response func- 
tion (RA-off-diagonal components) have been derived in 
Ref. [20. For the LDOS, we need the RA-diagonal com- 
ponents, <7f and g B , whose evolution is described by the 
coupled differential equations 



1 dg F 

2 dt 

1 dg B 

2 dt 



-Hg F 
-Hg B 



1 - A 2 dg 



Xf dXF 
A| - 1 %B 
dX B 



Ab 
Xf 



MAg, (B9a) 



\2 
£!_ 

Xf 



^MA 9 , (B9b) 
Ab 



where t is the distance from the point x measured in units 
of £, A.g = g B - g F , and 



M 



X F Xi 



1 



(Ab - Xf Y 



(BIO) 



in terms of two functions: 

g B X B ± g F X F 
= • 



(B13) 



Combining Eqs. (|B9a[) and (|B9b|) one immediately 
obtains^ that the function $_ exhibits a simple evo- 
lution with the Hamiltonian (|15|) : 



2#$_. 



(B14) 



Taking into account the boundary condition $_ = (Ab — 
X F )^(Q(x),s) at t = and the definition ([11]) of the 
Hamiltonian iJ, one arrives at 



$_(Q;M) - (A B -Af)*(Af,A B ;M), (B15) 



where 



*(A f ,Ab;M) =e- 2Ht *(g(.T),s) 



(B16) 



Evolution of the function $+(Q; i, s) is more complicated 
and can be in principle deduced from Eqs. (|B9p . 

E valuat ing the integral ([Blip using Eqs. (|B4[) . ([B5]) 
and (|B12|) . we get the LDOS for the symmetry class C: 

(p B (x)) =Rc$ + (l,l;t,s) + iR c y dA F J dX B 

P)p^—Sq{\b) 

x— *(\ F ,\ B ;t t s), (B17) 

oAb 

where the first term comes from the anomaly associ- 
ated with the singularity of the measure (|B4|) at the 
origin^ Though the general evolution of <£+ described 
by Eqs. (jB9|) is quite complicated, its value at the ori- 
gin, $_|_(A), can be readily found. Indeed, in the absence 
of coupling to the superconductor (Tj = 0), one should 
reproduce the normal-state result, (Q) — A. Hence 
$+(1,1; t,s) = 1, leading to Eq. ([2~T]|. 

Similar calculations for the symmetry class D [with 
Eq. |B6l) instead of (|B5[) ] leads to Eq. &o§. 
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Appendix C: Evolution at the Mott scale: 
and drift 



diffusion 



In this Appendix, we approximately calculate the wave 
function ^(Xp, Xb; t, s) defined in Eq. ([T31 for t ^> 1 and 
\n\ 2 <C 1. First consider the wave function 



*(A F ,A B ;0,s) 



-2Hs 



O 1. 



(CI) 



It interpolates between ^{Xf,Xb) = 1 at s = and 
*o(Af, A_b) given by Eq. (|T9|) at s oo. We will use the 
following properties of the wave function (|C1[) : 

(1) It does not develop any singularity at Xb e [l;oo) 
and X F G [-1; 1]. 



(2) It varies at the length scale A 



^> 1 in both 



Xp and As. In particular, the dependence on Xf is 
very weak. 

(3) It is universally normalized at the origin: 
^(1,1; 0, s) = 1. This follows from the supersym- 
metry (there is a zero mode of H with ^(1,1) = 1, 
and all the excited modes have ^(1, 1) = 0). 

Then we note that the fermionic part Hf of the Hamil- 
tonian (|16p has a discrete spectrum with a level spacing 
of order one, and therefore, at t ^> 1, only the lowest 
eigenstate of Hp survives. By the perturbation theory, 
its energy is of the order |/t| 4 , and, since we will be in- 
terested at the logarithmic scales t ~ ln|ft|~ 2 , it can be 
neglected, as well as its Xf dependence. Therefore, at 
t 1, we can approximate 



^(X F ,X B ;t, 



*(1,Ab;0,s). 



(C2) 



To understand the evolution in As, it is convenient to 
change to the logarithmic coordinate 9b defined by 



Xb = cosh 9 1 



(C3) 



In this variable, B ; 0, s) « 1 for 9 < 1ii|k|~ 2 (or 
even at larger 9b, for small s). In terms of 9b, the bosonic 
Hamiltonian takes the form 



4IH coth ^ + S cosh ^- (C4) 



This Hamiltonian consists of three terms. The first term 
describes a diffusion in 9b- The second term produces a 
drift towards the origin 9b = 0. Note that the velocity 
of this drift equals 1/2 at 9b 3> 1 and differs from 1/2 
only near the origin (in such a way that the total drift 
time changes by a small correction of order one). The 
third, potential, term of the Hamiltonian is relevant at 
9b > 6*m/2, where we define the Mott scale 



21n(| K | 



(C5) 



[we have introduced the coefficient 2 in the above defini- 
tion in order to relate 9m to the conventional definition 



Q of the Mott length scale Lm\. The third term sup- 
presses the wave function at 9b > 9 m /2 for t > 1, regard- 
less of the initial wave function (|C1|) . Therefore, at t ~ 1, 
the wave function ^(Af, Xb', t, s) may be approximated 
as a step function 



*(0b) « 0{0m/2 - e B ) , 



(C6) 



independently of Xf and s. Note that the width of this 
step function is of order one (in the units of Ob), but, at 
our level of precision, we are not interested in details at 
such a short length scale. 

Finally, we may approximately calculate the result of 
the evolution (|C2[) by applying the first two terms of the 
Hamiltonian (|C4[) to the step function (|C6|) . At our level 
of precision, we may also neglect the difference of the drift 
velocity from 1/2. We thus find the following result: 



^(X F ,X B :t,s) m - 



1 — erf 



9 B - 6 M /2 + t 



where 9b is related to As by Eq. (|C3 
In particular, at Xb — 1, we find 



2yft 



*(A F ,1;M) w - 



1-erf 



t - e M /2 

2Vt 



(C7) 



(C8) 



Note that the location of this erf- function is at Lm/2, 

which is twice closer than the similar "erf" feature in the 

correlation function of the local density of states in the 
wirej i4 i 50 

The result (|C8[) describes the crossover of &(Xf, 1; t, s) 
from one to zero around t w 9m /2 and is valid only in 
the vicinity of this point. To estimate the range of va- 
lidity of our approximations, we recall that we neglected 
the details of the evolution of the initial state at t < 1 . 
This effectively introduces an uncertainty At ~ 1 in the 
final result (|C8[) . Its range of validity may be therefore 
estimated from its insensitivity to this uncertainty, which 
leads to the condition \t — 9m /2| <C 9m (note that this 
range of validity is larger than the width of the crossover 
itself). 



Appendix D: Density of states at E 
class 



for the D 



The LDOS in the D class can be calculated exactly 
for quasiparticles right at the Fermi energy (E — > 0). 
In this limit, the "fermionic" variable q G [0, k] defined 
in Eq. (|2"0|) is small, and the wave function (fT3"|) can be 
expanded as 

*(p, q; t, s) = MP', t, s) + q 2 MPl t, s) + o(q 2 ). (Dl) 

Substituting this scries into the general expression 
([2"5| we see that the smallness of ^(1, As; t, s) — 
^(— 1, Xs;t, s) oc k 2 is compensated by the large value 
of the integral over Xb, leading to a finite limit at K — > 0: 



(p (x)) = l + 2 i/j 1 (p;t,s)pdp. 



(D2) 
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The function ^(p,q;t,s) is obtained through the 
Hamiltonian evolution (|1 3|) . Since the Hamiltonians H 
and H are related by the transformation (fl"6|). it suf- 
fices to study only the H evolution. In the order q 2 , the 
fermionic part (]17b[) of the Hamiltonian H acts as follows 
(at k = 0): 



o(, 2 



and one immediately obtains the evolution of an arbi- 
trary function xiPil) = Xo(p) + Q 2 Xi(p) + °(q 2 ) by the 
Hamiltonian H to the desired accuracy: 



-2Ht 



x(p, <l) = 1 + Q 



1-e 



-it 



8 

,2 -2ij 



UtXoip) 

+ q Z e- 2t U tX i(p)+o(q 2 ). (D3) 



Here 



U t = e 



-2H„t 



(D4) 



is the evolution operator for the bosonic Hamiltonian 
(fTTaj) in the limit k = 



B 8 \ p dp P dp P 



(D5) 



This Hamiltonian defined at the semiaxis p > has a 
continuous spectrum Ek = (fc 2 + l)/8 with the eigen- 
functions -Kjfc(p)/p. 

The general expression (flU)) for ^(p, q; t, s) can be iden- 
tically rewritten in the form 



•t>(p,q;t,s) = e- 2flt y (p,q) 
+ e- 2m {p 2 -q 2 )e 



^^Hs 1 - ^o(P: g) 

p2 _ g2 



, (D6) 



where 



*o(p) = *o(p,o), n (p) = a 9 2* (p,g)| 



F(p) 



1 - * (p) 



p 



Z(p) 



F(p)-n (p) 



p 



The evolution operator U t becomes trivial in the 
basis of the eigenfunctions of the Hamiltonian (|D5[) . 
where it amounts to multiplication by e~( fc + 1 )*/ 4 . 
Transition to the basis of Kik(p)/p is known as the 
Lebedev-Kontorovich transformation.— For the zero 
mode ^o(p,q) given by Eq. (|19p . wc obtain 



roo k 2 + I 

*o(p) = / tann 



irk K ik (p) 



P 



n (p)= / fed* — 

2 p kdk 



tanh 



7rfc -Kj fc (p) 



sinh 



2 P 
7rfc Kjk (p) 
2 p ' 



(D9a) 
(D9b) 
(D9c) 



Z{p) = - 



Wo (fc 2 + l)(fc 2 + 9) 



sinh^^M. (D9d) 



P 



In order to evaluate the action of the composite opera- 
tor U t p 2 U s in Eq. (|D8|) one has to reexpand pK^p) in 
the eigenfunctions of the Hamiltonian (|D5j) which is per- 
formed with the help of the Lebedev-Kontorovich trans- 
formation: 



pKik (v) = - Tidi ( /2 ' fc2 ) sinh7r; K «(p) (D10) 

4 2 Jo cosh 7rZ — cosh nk p 



Substituting Eqs. (|D9aj) and (|D9b|) into Eqs. (|D2|) and 
D7P , and integrating over p with the help of 



where we used the relation (fTS)) and extracted the zero 
mode of H . Now expanding in q 2 and successively ap- 
plying Eq. (|D3p . we arrive at 



^i(p;t,c») = 1 o e 2 V t ^o(p) + e- 2t U t U (p), (D7) 



,-2i 
1 _ e -2(t+s) 



U tP 2 U s Y(p) 



+ e- 2 ^U tP 2 U s Z(p) - e- 2t U t+s Y(p), (D8) 



/ K ik (p)dp = 
Jq 



2 cosh7rfc/2 ' 



(Dll) 



we obtain the result f4"2")l for the LDOS in the infinite 
wire, (p { °°\x)). 

The general expression for the zero-energy LDOS in 
a finite wire can be written with the help of Eqs. (|D2[) . 
(|D8|) . (|D9c|) . (jD9dl) and (|DT0|) as 



(po(x))=(pr ) (^))-2e 



tanh ^ e -( fc2 + 1 )( t + s )/ 4 
1 2 



oo />oo 



Jo 



fc 2 

kdkldl{k 2 



cosh7rfc 



I 2 ) . , . , ttZ 

— - sinh — sinh — 

cosh 7ri 2 2 



k 2 



e -2(t+«) 

fc 2 + 9 



-(fe 2 + l)s/4-(; 2 + l)t/4 



(D12) 
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